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In many high-profile applications, such as nuclear fusion and pumping of active media of short- 
wavelength lasers, it is necessary to achieve high specific input of power of an electromagnetic beam in 
a target. Diffraction sets the lower limit to the transverse dimensions of electromagnetic beams and 
represents a fundamental obstacle for electromagnetic heating of small or inaccessible regions. It was 
found, however, that it is possible to achieve efficient heating of cylindrical targets by electromagnetic 
beams with transverse dimensions that are several orders of magnitude greater than those of the 
cylinder. These counter-intuitive results have the following physical mechanism: the absorption in 
the cylinder causes a deep fall in the field distribution, and this fall causes diffractive diffusion of 
the field towards the axis from a large volume of the beam. The heating efficiency was rigorously 
calculated using the exact solution of the problem of diffraction on an infinite homogeneous cylinder 
(J.R. Wait, 1955). Several non-resonant domains of parameters were found that provide efficient 
absorption of the energy of a broad beam in a thin conducting cylinder. The typical asymptotic 
efficiency for very thin cylinders is (2-5) /L, where L=ln(w/d), where d is the diameter of the cylinder 
and w is the width of the beam waist in the longitudinal geometry or the wavelength in the transverse 
geometry. This efficiency is high even if the beam is several orders of magnitude broader than the 
cylinder. The relevant conditions are not too rigid and may be used to heat cylindrical targets to 
high temperatures in a number of important applications. External magnetic field may further relax 
the conditions of efficient heating. 

PACS numbers: 42.25.Fx;52.25.Os;52.38.Dx;52.50.Jm;52.50.Sw;52.80.Pi 



I. INTRODUCTION 

In many high-profile applications it is necessary to 
achieve high specific input of power of an electromag- 
netic beam in a target. No-contact methods are of spe- 
cial interest, as the resulting temperatures may be hardly 
compatible with any walls. Heating by electromagnetic 
waves is the most important of these methods. It is 
used for controlled fusion, for pumping of active media 
of short- wavelength lasers, in plasma chemical reactors, 
for material processing, and for many other applications. 
Diffraction sets lower limits to the transverse dimensions 
of the electromagnetic beam and thus presents a fun- 
damental obstacle for electromagnetic heating of small 
or distant regions. One may use shorter wavelengths to 
circumvent the diffraction constraints, but this is usu- 
ally connected with very serious technical difficulties and 
lower efficiency. Another possible way out is multiple- 
pass heating of small conducting regions - e.g., in res- 
onators or waveguides. But in this case heating is unac- 
ceptably slow for many applications, and it is practically 
impossible to heat distant regions (e.g., to form an ion- 
ized structure in the atmosphere for broadcast or radar). 
Non-linear effects, such as self-focusing, may be very im- 
portant, but they typically require an electron-density 
minimum in the region under heating and/or high beam 
intensity. 
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FIG. 1: Standard geometry. 1 - flat target; 2 - incident 
Gaussian electromagnetic beam. 



Thus, in standard geometries of target irradiation, only 
one dimension of the region under heating may be much 
smaller than the transverse dimensions of the electromag- 
netic beam (see Fig.^). On the other hand, it was found 
that in the longitudinal geometry (Fig. EJ) and the trans- 
verse geometry (Fig. [3}, if certain, not too rigid, condi- 
tions are met, it is possible to achieve efficient fast heat- 
ing of cylindrical targets by electromagnetic beams with 
transverse dimensions that are several orders of magni- 
tude greater than those of the cylinder (Ref. Q, Q, 
0, |g, H, 0, 0, @). These counter-intuitive results 
have the following physical mechanism: the absorption 
in the cylinder causes a deep fall in the field distribu- 
tion, and this fall causes diffractive diffusion of the field 
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FIG. 2: Longitudinal geometry (the axes of a Gaussian beam 
and a cylinder coincide; the energy flows mainly along the 
axis), f - cylindrical target; 2 - incident Gaussian electro- 
magnetic beam (energy flows from left to right). 




FIG. 3: Transverse geometry (the axes of a converging cylin- 
drical wave and a cylinder coincide; there is no energy flow 
along the axis; the wavelength is considered as the transverse 
dimension of the beam). 1 - cylindrical target; 2 - incident 
converging cylindrical electromagnetic wave . 



towards the axis from a large volume of the beam. The 
effect is linear (for given properties of the cylinder, the 
heating efficiency does not depend on the beam power) 
and scalable (may find applications in various frequency 
ranges). As these results are rather surprising, much 
effort was made to reliably confirm them. For instance, 
for the longitudinal geometry, three methods were used: 
a qualitative (" back-of-the envelope") estimate, a semi- 
quantitative ("one- wave") calculation, and a quantitative 
calculation, where the Maxwell equations were solved vir- 
tually without any approximations. The latter was used 
to confirm the results of the semiquantitative approach, 
which was used to derive an asymptotic formula for heat- 
ing efficiency and determine the parameter domains that 
provide high efficiency. 

In the qualitative estimate (Section llVI) . it is assumed 
that the values of the frequency of the electromagnetic 
field and the conductivity of the cylinder are such that 
the skin-depth and the diameter of the cylinder are of 
the same order of magnitude. Then it is assumed that in 



this case the refracted electric field inside the cylinder has 
the same magnitude as the electric field in the incident 
beam (the electromagnetic field penetrates the cylinder 
to the skin-depth). Then we obtain that the Joule heating 
power is of the same order of magnitude as the power in 
the incident beam. 

Of course, diffraction problems are notorious for strong 
dependence on the parameters, and this primitive esti- 
mate is often incorrect, but it gives reasonable order-of- 
magnitude results in some important cases and demon- 
strates the scaling of the problem. For example, if we 
reduce the diameter of the cylinder diameter by half, 
the optimal conductivity increases fourfold (as the skin- 
depth is inversely proportional to the square root of the 
conductivity), and the absorbed power remains roughly 
the same. Alternatively, if the width of the beam waist 
increases twofold, the electric field in the cylinder falls 
foufold, but the length of the beam waist increases four- 
fold, so the absorption efficiency does not change. 

In the semiquantitative approach, it is assumed that 
the refracted electric field inside the overlap of the cylin- 
der and of the Gaussian beam waist is the same as the 
refracted electric field for a cylindrical wave that has the 
field intensity on the axis equal to that in the center of 
the Gaussian beam and the width of the first maximum 
equal to the width of the beam waist. Then the exact so- 
lution of the problem of diffraction of a cylindrical wave 
on an infinite homogeneous cylinder is used (Ref. [Tol]). 

The quantitative calculation was performed as follows: 

1. For rigorous calculation of diffraction of a Gaus- 
sian beam on an infinite homogeneous cylinder, exact 
solutions of the free Maxwell equations were found (Sec- 
tion^ that are approximated by Gaussian beams with 
great accuracy when the beam waist radius is much 
greater than the wavelength (if this condition is not ful- 
filled, a Gaussian beam is not a good approximate solu- 
tion of the free Maxwell equations). This result may be 
important for many applications. 

2. For the cylindrical waves, the Wait's solution was 
calculated (in this exact solution of the Maxwell equa- 
tions, each cylindrical wave incident on an infinite ho- 
mogeneous cylinder generates two reflected and two re- 
fracted cylindrical waves). 

3. The resulting refracted fields were calculated by 
integration of the fields of the refracted cylindrical waves 
using fast Hankel transform. 

4. The heating efficiency was calculated as a ratio of 
the power of the Joule heating in the cylinder (within the 
length of the waist of the Gaussian beam) to the power 
in the incident free Gaussian beam. 

Thus, virtually no approximations were used in this 
rigorous calculation. The results of the semiquantitative 
approach coincided with those of the rigorous calculation 
with 50% accuracy. 

Let us describe the structure of this first part of the 
work. 

In Sectional exact solutions of the free Maxwell equa- 
tions are described that are approximated by Gaussian 
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beams with great accuracy when the beam waist radius 
is much greater than the wavelength, and their asymp- 
totic properties are proved. It is also proven that these 
solutions are normalizable. 

In Section IIIII the exact solution of the problem of 
diffraction of a cylindrical wave on an infinite homoge- 
neous cylinder (Ref. |lOj |) is summarized. 

In Sect ion llVl the qualitative and the semiquantitative 
approaches are described. In the framework of the lat- 
ter, an asymptotic formula for heating efficiency in the 
longitudinal geometry is derived. 

In Section an asymptotic formula for heating effi- 
ciency in the transverse geometry is derived. This cal- 
culation is relatively simple and easy to verify, but has 
many common features with the more complex calcula- 
tions for the longitudinal geometry. 

In Section EH some conclusions are formulated. 

Thus, this first part contains considerable technical 
material. In the second part, the parameter domains are 
determined that provide efficient heating in the longitu- 
dinal geometry. The relevant conditions were formulated 
in part in Ref. JUi 0! they refine and generalize the 
results of Ref. The conditions are confirmed by rig- 
orous calculation, and some applications are discussed in 
the second part. 



II. EXACT SOLUTIONS OF THE FREE 
MAXWELL EQUATIONS AND GAUSSIAN 
BEAMS 

It is highly desirable to validate the results of the 
present paper through comparison with exact solutions 
of the Maxwell equations. However, Gaussian beams are 
just approximations to exact solutions of free Maxwell 
equations. So in this section the author informally de- 
rives some exact solutions of the free Maxwell equations 
and proves that Gaussian beams approximate them with 
high accuracy. These solutions were found in Ref. @- 
In Ref. 01 > drawbacks of the earlier approaches to this 
problem are discussed, and other beam-like solutions of 
the free Maxwell equations are found. 

We consider a monochromatic beam with frequency 
lj and use such a system of units that the wave vector 
in free space fco = — = 1. Solutions of the Maxwell 
equations in homogeneous media with permittivity s and 
magnetic permeability [i are described by the electric and 
magnetic Hertz vectors II and II'. The Hertz vectors 
satisfy equations: 



An + /c 2 n = o, 



(1) 



e = - v x (v x n) + ik v x n'. (4) 

e 

In the cylindrical coordinates {/?, tp, z} or in the cartesian 
coordinates {x,y, z} we can choose the Hertz vectors in 
the following form: 



n = {o,o,n}, 



IT = {0,0,11'}. 



(5) 



(6) 



In free space e = fx = 1. To emulate a Gaussian electro- 
magnetic beam with circular polarization let us assume 
that in the symmetry plane of the beam (z = 0) we have 



E(x,y,0) ={z,-l, 0}exp(- 



2S 2 



H-(x,y,0)={M,0}exp(- 



x 2 + y 2 , 
26 2 ■ 



(7) 



(8) 



It should be noted that, for the sake of simplicity, some 
conventions are used here, so, the actual electric and 
magnetic fields are equal to the real parts of the prod- 
ucts of the expressions in Eqs. Q and © by the time- 
dependent factor exp(— iut). Let us write the Fourier 
expansions of Eqs. Q and JHJ : 



S 2 f ( ^kl+k 2 , 
E(x, y, 0) = {i, -1, 0}— J cxp ( -5 2 ^-JL 

x exp(i(k x X + k y y)) dk x dk y , 



f) 2 r ( k 2 + k 2 
H(x, y, 0) = {1, i, / exp -5 2 ^-^ 



(9) 



x exp(i(k x x + k y y)) dk x dk y 



(10) 



However, these expressions are not compatible with 
the free Maxwell equations, as, e.g., on the one hand, 
H z (x,y,0) = 0, therefore H z (x,y, 0) = and (V X 
E) z (x,y,0) = 0, and the latter is not true. So we have 
to modify these expressions. The Fourier expansions sug- 
gest that the fields may be constructed as a superposition 
of plane waves. Let us consider a plane electromagnetic 
wave with a wave vector {0, 0, 1} - it propagates in the 
direction z. To obtain a plane wave with a wave vector 
k x ,k y ,k z , where k 2 + k 2 + k z = 1, we may multiply the 
field vectors by the matrix 



An' + k 2 u' 



0. 



(2) 



where k 2 = . The electric and magnetic fields E and 
H can be calculated from the Hertz vectors as follows: 

H = -4Vxn + -Vx(Vxn'), (3) 

A 4 




l + k z 



k k\ 






Ky I 


+- 




■y k z J 






k 2 


k x k y 





ky 


k x 















(11) 
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and change the phase factor. Following such reasoning, or, if we change the variables p — > A, r — > p, a — > ip, 
we obtain (or, rather, choose) 

n (p, tp,0) = w%exp(i<p) [ exp ( -w^ )Ji(Xp)dX 



E z (x,y,0) = ^ J exp ( -to 2 



x exp(i(fc x x + k y y)) dk x dk y , (12) 



It can be shown that 

9 
10, 



H(p, if, 0) = — exp(i^) 1 - exp - 



(19) 



(20) 



H z (x,y,0) = -iE g (x,y,0). 
As E z = -A 2 n = -{d 2 x + sg)n, 



(13) 



We need to know the fields for all values of z. We 
choose the following expressions for z-components of the 
Hertz vectors: 



n(z,y,0) 



2tt 



/ 2 fc 2 + fe 9 
exp -w Q — — 



+ ky 

1,2 I 1,2 



n(w) = 5 



ex P - w 7T- 



1 



ik x -\- k u 



x exp(i(fc x x + dk x dk y 



2tt 



( 2 + ^ u 

ex P -Wq JT - - 



1 



x exp(i(k x x + k y y)) dk x dk y . (14) 
Let us change the variables: k x = pcosip, k y — psimp: 

,2>r 



W = -ill. 



x exp (i (k x x + k y y + yd - fc 2 - k 2 zj) 

x 

(21) 



U(x,y,0)= I exp [ -w^ 



2 / p(sin + z cos </?) 
x exp(ip(x cos (y3 + ysmip))dpd(p = 

2\ ! 



— I exp I -w5 — 



27T 



^z 9 



We have limited the region of integration by the condition 
k x + ky < 1 to eliminate reactive fields and leave running 
waves only. In Eq. (|21|l the solution is presented as a 
linear combination of plane waves. For our purpose it is 
preferable to present the solution as a linear combination 
of cylindrical waves: 



2 J iexp(—itp) 



x exp(?p(x cost/? + y sin (p))dpdip. (15) 

In the cylindrical coordinates (x = rcosa, y — rsina), 
as cos a cos ip + sin a sin ip = cos(a — <p), we have: 



f 1 ( A 2 
n (p, ip, z) = Wq exp(i<p) J exp ( -w% — 

(iy/l - A 2 z) Ji(Ap) dA, 



x exp 



n' = -m. (22) 



n(r, a, 0) = — / exp ( -w 2 — J / (-i exp(iy)) Let us find the asymptotics of Eq. J22J| for the case w > 
., 1 (when the Gaussian beam waist width is much greater 

X exp(ipr cos(a - ip)) dpd(p. (16) tha.n the wavelenethV Let, us assume tha,t a is a. laxem 



x exp(ipr cos(a — <p)) dpd<p. 
Let us calculate the following integral: 



exp(ipr cos(a — ¥>))(— iexp(i<p))d<p 



i / exp(ipr cos(a — ip)) e~xjp(i((p — a)) exa{ia)dip 
= — iexp(ia) I exp(ipr cos ip) exp(i<p)d(p 



(when 

than the wavelength). Let us assume that a is a large 
parameter . In view of the scaling of Gaussian beams (the 
beam waist length increases as the square of its width), 
we assume that too = aw' , p — ap' , z — a 2 z', where w' Q , 
p', and z 1 are fixed, and a — > oo. 



H(p,<p,z) = o^w'q exp(i<p) I exp ( — Wq 2 q; 2 



A 2 



x exp (iy/l - A 2 a 2 z'^j Ji(Xap') dX. (23) 



= — iexp(ia) / exp(iprcos^) exp(itp)d<p = 
Jo 

= — i exp(ia)27riJi(pr) = 27rexp(ia) J\{pr) (17) 
(Jfe(x) is the Bessel function). Thus, 



II(r, a,0) = u)Qexp(ia) / exp ( — Wq^— I J\(pr)dp 



A 2 
2 



(18) 



i-^ + VT - A 2 " 



2 ^l-f + x/I - ^ 



(24) 
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therefore, 



H(p, <p, z) = a w exp(itp) / exp — w a 



■ cxp ( i [ 1 — ) <rz' 



x exp 



■a 2 z' Ji{Xap) dA.(25) 



4(1-^ + ^1^) 

Let us change the integration variable: x = Xap' . 

II(p, ip, z) — a 2 w' 2 exp(j <p) 

,2 / 



exp(ia z ) exp — 



2/ 



exp 



/2 



x exp 



A —2 / — 4- / 

x a p z 



4 1 



2p 



ap' 



-2- exp(i(^) exp(ia 2 z') 



x exp — i- 



rup 

x / exp 
Jo 



A — 2 /— A t 

x^a z p ^z 



4 1 



2a 2 p'- 



Ji(x)dx. (26) 



Let us assume that p' ^ and a 1 / 4 < ap' (a is great 
enough). Let us also note that 

A 2 



1/2 < /(A) = 1 - — + \/l-A 2 < 2, (27) 

on the segment < A < 1, as /(A) is monotone decreas- 
ing on this segment, and its values at the ends of the 
segment are 2 and 1/2. We also note that | Ji(x)\ < 1 for 
real x. 

The last integral in Eq. (|26|l may be presented in the 
following form: 



ap 

.1/4 



(28) 



We have: 



a p 

,1/4 



< 



< 



exp 



,1/2 



2p> 



-K 2 + ^) 



ap — a 



1/4 



(29) 



for 



x(l + 0(a" 1 ))Ji(x)dx, (30) 



,1/4 



/ ex p(-^72K 2 +^')) xO(«- 1 )Ji(x)dx < 

< lO^-^la 1 / 4 ^0 (31) 



for 



exp 



2/ 



(u^ 2 + iz') ) x Ji(x) dx = / +o(l), 



exp 



at ITS ■ i\ J l( x ) dx 
2p' z (w + iz ) / 



1 — exp 



2(w' Q 2 + iz') 



(32) 



(cf. Ref. [jj, P-186, formula 2.12.9.2). Therefore, for 
large a 
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hup, tp, z) w exp(i(p) exp(za z ) 

P' 

x(l-exp(- 2( ^ + . z/) 



= — exp(z r 9) exp(iz) ( 1 — exp ( — 



2(w 2 + iz) 



(33) 



We also have 



n' = -m. 



In free space e = 1 and = 1, therefore, 



(34) 



H = -iV x U iV x (V x II), (35) 



# = v x (v x n) + v x n, 



thus, 



E = iH. 



(36) 



(37) 



Taking into consideration Eq. (J5J), we have in the cylin- 
drical coordinates: 



_ , 1 an an n , 

Vxn= 



(38) 



6 



vx(Vxn) 



dz \ dp J ' dz \p dp J ' p \dp \ \ dp J J dp \p dp J 

r 



(39) 



an =m 

dp 

Let us present II in the following form: 

n = ni + n 2 , 

where 



(40) 



(41) 



dlli 

dz 



ini, 



(45) 



III = — cxp(z</?) exp(zz), 



(42) 



dll 1 _ 1 
5p P 



(46) 



w 2 ( 
n 2 = - exp(itp) exp(iz)exp ( - 

P V 



2(wg + iz) 



We have: 



dIL_ 



= *ni, 



• (43) 



(44) 



vxn 1 = | -n 1 ,-n 1 ,o;>, 
p p 



(47) 



v x (v x ni) 



dz (p Hl ) ' 9z (p Hl ) ' p Gp 9< p (p Hl 

i 11/1 1 

-III, — III,- — III + -III 

p p p\ p p 



(48) 



therefore, LTi does not contribute to the fields. 
dU 2 



dtp 



= *n 2 , 



(49) 



dU 2 -p 

— — - — llo 

dp 



Wq + iz p \Wq+iz p 



-n. 



i 



2 



Wq +IZ p 



,(51) 



an 2 

dz 



= m 2 + n 2 



p 



= *n 2 i + 



2(w 2 +iz) 2 

„2 



2( w 2 + i Z )2 y ' 



(50) 



v x n 2 = <j -n 2 , . 



p 



i 



WO + p 



- n 2 ,o , (52) 



(v x (v x n 2 )) p = — ( -n 2 j 



p - + - ) ) = -ffl 2 ( i + 



wg +iz p 



= -iIL 2 



2 ~. ■ ■ ~ ) II .' ~ ~. = 



2(w 2 + iz) 2 J \WQ+iz p) (w 2 +iz) 



wl +iz p 2(w 2 l + iz) 3 2{wl + iz) 2 (w 2 , + izf 



, (53) 



(v x (v x n 2 )) v = |- f-n 2 ) = -zn 2 ( 1 + / , 

<9z VP / P V 2 ( w o + 2Z ) 



= -112 



n 2 - 



p 2(w 2 + i z )- 



(54) 
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(V x <v x n 2 », . 1 (| (n. ( ' 



-n. 



1 



p 



i 



p 2(w 2 +iz) 2 / p 

2 

n 2 ' 



n 2 



p 



i 



l 

— + - 

iz p 

1 



(wq + iz) 2 Wq + iz Wq + 



P 



Wq + iZ 



1 



n 2 



-l 

2p 



i 



Wq + iZ p 



-n 2 



ZZ Wq + IZ P Z 



n 2 1 = 

n 2 

P 2 



(w 2 + iz) 2 



(55) 



iU 2 



P 



w' 2 + iz 2(wq + iz) 3 2(w 2 + iz) 2 (w 2 

2 \ 1 



,n 



p 



= — w exp(itp) exp(iz) exp 



P 



2(u>q + iz) / Wq + iz 



w 2 + iz 2(wq + iz) 2 
2 1 

1 - 



v x n 2 + v x (v x n 2 

„ 2 

,n 2 



p 



i 



(w 2 , + iz) 2 
P 



2(w 2 + iz) 2 2(w 2 + iz) / 2(w 2 +iz)' w 2 + iz 



• (56) 



If (f — 0, z = 0, p — 0, wo 3> 1, the fields have the 
following components: 

£ = VxII 2 + Vx (V X II 2 ) = {i,-l,0}, (57) 



are valid for the Hertz vectors of Eq. 1221) . From Eqs. 1221) 
and H24|l we have: 



n(p, <p, z) = wl exp(i<p) / f(X)J 1 (Xp)dX, (61) 
Jo 



H = -iE = {l,i,0} 



(58) 



With an accuracy of o(l) the field components may be 
presented in the following form: 



E 



x exp 



w, 



^— exp (i ip) exp(iz) 



■ iz 

,2 



2(w 2 + iz) 



0,-1,0}, 



(59) 



H 



x exp 



w. 



o 



Wq + iZ 



2(w 2 +iz) 



exp(itp) exp(iz) 



{l,i,0}, 



(60) 



(The field components are not small in comparison with 
1 only if p < \J\wq + iz\; we also take into account that 
wo » 1). Thus, the fields have the space distribution 
typical for a Gaussian beam. However, so far the fields 
were calculated from the Hertz vectors of Eqs. <|33[) and 
(|34|l . so let us prove that the same asymptotic expressions 



where 



f(X) = exp 



x exp iz 1 



2 4 (l - f + x/T^A 2 " 




We need to calculate the derivatives of n(p, <p, z 



an 

dip 



(62) 



(63) 



and to obtain §2, |P an d §4r we need to make the 

£72 ' op ' ap z 

following replacements in the right-hand side of Eq. I|61|l : 



/ \2 

/(A)-i 1-^-- 



A 4 



4 (l - 4j- + Vl^A 2 ' 



/(A), 



Ji(Ap)->AJ{(Ap), 
Ji(Ap)-> A 2 Jf(Ap),(64) 



correspondingly. We have: 



.e = v x n + v x (v x n) = ^ -— + 

1 p o<p 



tan a 2 n 1 a 2 n 



dzdp' p dzdip 
.P 



an i 

dp' p 

i an 

dzdp 1 p dz 



an 

~ap" 
an 

ap ' 



a 2 n i a 2 n 



^ dp 2 

ian 

p ^p 



pdip 2 

a 2 n 



ap 2 



P 



(65) 



r 



The terms with n in the last line of Eq. I|65|) may be esti- 



mated in the same way as n in Eq. (|23|l and further. For 
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the terms with and we should make the following 
replacement: 

Ji(Ap) - {o,-AJ{(A P ),~J((A P ) - A 2 J('((Ap)| . (66) 
In Eq. H26|) we make the replacement: 



Ji(x)-> U~4(x),~Ji(x) 
{ ap' ap' 



a 2 p' 2 



Jf (*) . (67) 



The estimates closely follow those in Eq. (|26ll and further. 
The remaining integral analogous to that in Eq. 13 2|) 
may be calculated by differentiation of the latter with 
respect to a parameter. For the terms with and 
in Eq. I|65ll we make the replacement of Eq. (|64|l for /(A) 
and the following replacement: 



J 1 (Ap)^AJ((Ap),~J 1 (Ap),0 
P 



(68) 



To estimate these terms we make the following replace- 
ments in Eq. I|26|) : the factor 



9 = i\ 1 



4 —4 —4 

x a p 



2a 2 p n 



4 1 



2a' 1 p' 



(69) 



is added, and 



I ap' ap' J 



(70) 



Then we follow the estimates of Eq. I|26|) and further. 
The remaining integral has the following form: 



1/4 



exp 



x <^ — jJ'i(x), 7 Ji(x),0> dx. 

1 ap' ap' 1 



(71) 



The first two terms in the expression for g (Eq. JSHJ) give 
the contributions that are estimated in the same way as 
the integral in Eq. (1.'S2[| and present its derivatives. For 
the remaining integral 



/ = 



y l/4 



exp 



x 
2~p 



n [w\ + iz'\ 



4 _4 _4 

x a p 



4 1 



2a 2 p' ; 



jJ'^^Mx)^} dx (72) 
ap' ap 1 



we have |/| < a _3 o(l), so it does not make significant 
contributions to the derivatives. 

Now let us prove that the solution of the free Maxwell 
equations (Eq. is normalizable, i.e. integrals of 

products of any two field components (such as E^Hp, E 2 
etc.) over any transverse section z = const are finite, so 
energy flow through a transverse section and total elec- 



tromagnetic energy between two transverse sections are 
finite. It may be seen from Eq. I|22|l that all field compo- 
nents are continuous, so it is sufficient to prove that for 
any z field components are 0{p~?) for large p. Eq. I|37|l 
implies that it is sufficient to prove this for components 
of electric field. We obtain from Eq. I|65|l : 



E{p,(p,z) = w 2 J exp (i(p + i\J\ - A 2 z - w 2 y^j 

x (-Ji(Ap) +iy/l- X 2 XJ[{Xp),-XJ' 1 (Xp) - - Vl - A 2 Ji(Ap), --AJ((Ap) - X 2 J{'(Xp) + \ji(Xp) } d\. (73) 
[p p p p 2 J 



If we denote exp ( i(p + - X 2 z - wf t \- J as q, then Therefore, the term with p is 0{p ) for large p. Let 

us now assess terms with p~ l , using the first term as an 

M<i, 
\MXp)\ < 1, 

\MXp)\ < 1. (74) 



9 



example: 



qJi(Xp) dX = 



--qJo(Xp) 
P 



i ! ,1 



— / qJo{Xp) 
P Jo 



o + PJo 

~iXz 



-d(J (Xp)) = 
P 

J Q (Xp)dq^O(p- 1 ) 
w^X) dX = Oip' 1 ) (75) 



(the singularity in the last integral is integrable). The 
other two terms with are assessed in a similar way. 



Let us assess another term: 



P Jo 



Ji{Xp)q 



qVl-X 2 XJ[{X P ) dX = / qy/l-X*X 

I! JO 

iXz ,,\ r rrr . A 2 



d(Ji(Xp)) 



((^-^)yr^__£_ + % /r^. 



(76) 



The last integral may be presented as a sum of two inte- 
grals: 




The first integral may be assessed using another integra- 
tion by parts, similarly to the previous integrals (this in- 
tegral does not contain any singularities). The second in- 
tegral is 0(p~2 ) ; as for A > i we have Ji(Xp) — 0(p~ 5 ). 
Thus, the whole term is 0(p~^). 

The remaining two terms can be assessed in a similar 
way. 



III. THE EXACT SOLUTION FOR THE 
PROBLEM OF DIFFRACTION ON AN INFINITE 
HOMOGENEOUS CYLINDER 

We shall extensively use the exact solution of the prob- 
lem of diffraction on an infinite homogeneous cylinder by 
J. Wait Ref. 0], so ^ us outline it here. Let us consider 
cylindrical waves of TM type: 

n = Z n (Xp) exp(i(jz + mp)), 
H = l^^-Z n (Xp),ik o XZ' n (Xp),0^ exp(i(7z + nip)), 

E=\ l -^Z' n {Xp),- 1 —Z n {Xp),-Z n (Xp)\ 
ye ep e J 

x exp(i(-fz + nip)), (78) 



and TE type: 

If' = Z n (Xp) exp(i(7_z + nip)), 

H' = { l J± Z ' n (Xp),-^Z n (Xp), -Z n (Xp)\ 
I P PP P J 

x exp(i(7z + nip)), 
E' = ^-^yZ n (Xp),-ik XZ' n {Xp),ojexp(i{ 1 z + nip)), (79) 

where Z n is a cylindrical function (J„ or H^), and 

A 2 = fc 2 — 7 2 . (80) 

The cylinder is determined by the equation p < a, so a is 
the radius of the cylinder, the permittivity and the mag- 
netic permeability are E\ and p\ outside the cylinder and 
Ei and p2 inside the cylinder. Let the incident field be a 
linear combination of cylindrical TM and TE waves with 
the following z-components of the electric and magnetic 
Hertz vectors, correspondingly: 

uo = aJ n (Xip)F n , 

v =pj n (X lP )F n , (81) 

where F n = exp(i(7z + nip)). Then the reflected field 
is also a linear combination of cylindrical TM and TE 
waves: 

Mi = a 1 ^ l H^ 1 \X 1 p)F n , 
v 1 = b l , n H^(X 1 p)F n , (82) 

and the same is true for the refracted field: 

U2 = a 2 , n Jn(X2p)F n , 

V2 = b 2 ,nJn{X2p)F n , (83) 

(ui and u 2 are the z-components of the electric Hertz 
vector, v\ and v 2 are the z-components of the magnetic 
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Hertz vector). The coefficients a\ >n , a2,m ^i,ni b 2 . n are 
evaluated from the continuity conditions for the tangen- 
tial components of the fields on the surface of the cylin- 
der. We shall omit subscripts n and superscripts (1) for 
simplicity. For p = a we have (superscripts inc, rfl, and 
rfr relate to the incident, reflected, and refracted fields; 
Aia = pi, X 2 a = p 2 ): 



Tjtnc 



i 772 , 

aik \iJ'(pi) — P J{pi) + aiikoXiH'ipi) — 

api 

-bx^Hfa) = ik X 2 a 2 J'{p 2 ) - ^b 2 J(p 2 ), (84) 



api 



a[i 2 



7n 
E\a 



aJ(pi) - ifc Ai/?J'(pi) 



771 / El X\ J(p 2 ) J( Pl ) 

' e ia \e 2 xlH{pi) a2 n{ P x) a j 



) H( Pl ) 



., , f Mi Aj J(p 2 ) Jjpi) A ,. 
-tk Xi —^ jj, M - J 3 H (pi) = 
V/i2AfiJ(pi) H{ Pl ) J 

= -^a 2 J(p 2 ) ~ ik X 2 b 2 J'(p 2 ). (91) 

e 2 a 

Thus, we have a system of two linear equations for a 2 
and b 2 : 



c\\a 2 + c 12 b 2 = ci, 



(92) 



aJ(p 1 ) - ik XiPJ'(pi) a 1 H(p 1 ) - 



-ik XibiH' (pi 



771 
£20 



a 2 J{p 2 ) - ik Q X 2 b 2 J'(p 2 ), (85) where 



HI 



fjrfl _ jrrfr 



Ml 



/?J(pi 



A U oY„ N _ A 2 



Mi 



-hH( Pl ) 



M2 



&2</(P2), 



We obtain: 



61 



Mi 



X{H{p 1 ) VM2 



A '&2J(P2) - ^J(pi! 



Mi 



Mi A| J(p 2 ) , J(pi) 



fi 2 X(H( Pl ) Hipx) 



a 1 



£ i ( ^2 t / \ Af 



— o 2 J(p 2 ) - — aJ{px) = 



£1 



£l A| J(p 2 ) J(px) 



e 2 X\H{ Pl ) a2 H( Pi y 
Let us exclude a\ and 61 from Eqs. (|84I85|) : 

aik Xi J'{pi) - P J(pi) 



+ 



£1 A| J{p 2 ) J(pi) 



-a 2 



a I ikoXxH'ipi) 



s 2 X\H(p i y H( Pl 
Mi A| J(p 2 ) , </(Pi) A 7" 



\x 2 X{ H(pi) H(pi) J api 



P ^H( P1 



(86) 



M a j(p x ) + ^aiHfa) = ^-a 2 J( P2 ). (87) 
£1 £1 £2 



/?, (88) 



(89) 



ik X 2 a 2 J'(p 2 ) - ^-b 2 J(p 2 ), (90) 
a/j 2 



c 2 ia 2 + c 22 6 2 = c 2 , (93) 



1 1 e 2 J'{p 2 )\ £ip 2 p 2 , . 

ui \ 77 — r ik — J(p 2 ), 

pi H(pi) p 2 eiJ(p 2 )J £ 2 a 



C12 



pt p 2 



a/i 2 



C21 



C22 



Pi P 2 ) ae 2 
1 fTfa) 1 ^2 J'fe) 



pi H(pi) p 2 pi J(p 2 

v (-M1P2) ., P2 w , 

x ik — J{p 2 ), 

M2 a 

/ g^i) J'(p x )\ 

Cl = m " oAl l^ij- T^yJ J(pi) ^ 

aifcoAi 2i 

-^bl); 



c 2 = ik XiPJ(pi) 



H(p 1 )J(pi) npi 
fJ'( Pl ) H'( Pl ) 
\ J(pi) H(px) 
ikoXiPJ(pi) ( — 2i 



J{pi)H(pi) \irpx 



(94) 



We have used the well-known formula for the Wronskian 
of cylindrical functions: 

W[J n (z),H n (z)] = —. (95) 

Let us solve the system using the Cramer's Rule. 

£1 ik plJ(p 2 ) fix ik Q p^J(p 2 ) 



D' = 
where 



cn c i2 
c 2 i c 22 



£2 



M2 



-D, (96) 
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D 



1 g'(Pi) 

Pi H(pi) 



1 £2 



J'(P2 



_yn_ ( 1 

ifc £i 



7" [ 1 

ikoHi \p( 



1 g'(Pi) 



1 

U2 J' 



P2 £1 J(P2) 

pj) I /'! .'/(/'i i ;<j /'i •/!,._.; 

1 1 £ 2 J'(p 2 )\ / 1 ff'^l) 



9 2 

7 n 



2 

Pi 



1 M2 J 7 (pa) 

Pi H{pi) p 2 e 1 J(p 2 ) J Vpi h (Pi) P2M1 Jfe) 



(97) 



Then we have: 



£>i 



Cl C12 
C2 C 2 2 



__ 2a fcp 

Q7T H{pi) 

2/3 fcp 

Q7T -ff(pi) 



7" [ 1 

*0Mi \pf 



1 g'(Pi) 
pi ff(pi) 



1 

1 M2 J'(P2) 
P2 Ml ^(P2) 



Mi ., v\J(vi) 

— ifc 

M2 a 



7" / 1 
ikofJ.i \pf 



1 g'(pi) 
Pi -ff(pi) 



1 

1 M2 J'(P2) 
P2 Ml ^(P2) 



2k 



airH(pi) \x 2 



Mi v\J{V2) 
ik . 



hence, 



02 



A. 



-2/: 



£2 



nH(pi)D e\ p 2 J{p 2 



a- 



£2 1 



-1% 



1 ^'(Pi) 



1 H'fpt) 1 M2 J'(p 2 ) 



Pi H(pi 

1 /Lt 2 J'(P2 



e\ J{p 2 ) itp 2 H(j>\)D \Pi -ff(pi) P2 Mi ^(P2) 



M2 Ml J{P2) 

£2 1 



■0 



fcoMi 
2 771 



1 

p| 



£1 J(p 2 ) irpiH(pi)D ^oMi 



Pi 



2 

Pi 



9 
T>2 



(99) 



Similarly, 



#2 = 



Cn Cl 

C21 c 2 



1 g'(pl) 

Pi f"(pl) 



x 

ik, 



oei \p{ 



1 E 2 J'(P2) 

P2 ei J(p2) 

4 



ei ik plJ(p 2 ) 2fc 
e 2 a anH(pi) ' 



(100) 



L> 2 -2i /i 2 1 { ijn { 1 1 \ / 1 F'Oi) 1 e 2 J'(p 2 



£>' nH(pi)D ^ipiJ(p 2 ) V fc o£i VPi P2/ VPi #(pi) P2 £1 >/(P2 



g M2 1 2 7 n / 1_ _ 1_\ | Mg 1 -27 / 1 g ; (pi) _ 1 £2 J'{P2) 

Vi J(P2) irplH(pi)D k ei \j% pj) m J{p 2 ) nplH( Pl )D \ Pl H( Pl ) p 2 e x J(p 2 ) 



IV. THE CONDITIONS OF EFFICIENT 
ABSORPTION OF A BROAD GAUSSIAN 
ELECTROMAGNETIC BEAM IN A THIN 
CONDUCTING CYLINDER 

Diffraction problems are often highly sensitive to vari- 
ations of their parameters. Therefore, confirmation of 
the results of any calculations by rigorous solution of 
Maxwell equations is highly desirable, and this proce- 
dure is discussed in the next section. However, one also 
needs simpler approximate procedures to make estimates 
and to determine the most promising parameter domains. 
Therefore, in this section the author describes a sim- 
ple (and sometimes inaccurate) qualitative estimate of 



the part of the power of a Gaussian beam absorbed in 
a conducting cylinder and a semiquantitative one-wave 
method. The latter is used to deduce an asymptotic for- 
mula for efficiency of the absorption and to determine 
the parameter domains where efficient absorption takes 
place. The results of the semiquantitative methods were 
confirmed by rigorous solution of the Maxwell equation. 
This section mostly deals with the longitudinal geome- 
try, as for this geometry conditions of efficient absorp- 
tion of a broad beam in a thin cylinder were not known, 
so the author investigates this case in much detail. For 
the transverse geometry some resonant cases of efficient 
absorption were known (see, eg.,Ref. [13), so the author 
describes just one non-resonant parameter domain where 
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efficient absorption takes place in a cylinder with a ra- 
dius that is much less than the wavelength. This result is 
important as it complements the results for the longitudi- 
nal geometry. To describe the qualitative estimate for the 
longitudinal geometry, let us consider a Gaussian beam 
of electromagnetic waves with frequency lu in free space. 
Let the length of the conducting cylinder be greater than 
the beam waist length l± denned by twofold field inten- 
sity attenuation, its radius a be (much) less than that of 
the beam waist r± defined by e-fold field intensity attenu- 
ation, and its conductivity a be such that the skin-depth 
be roughly equal to the cylinder diameter: 



4a 2 = 



2irui<j 



(102) 



Let us place the cylinder in the beam in such a way that 
the axes of the cylinder and the beam coincide, and the 
cylinder overlap the whole beam waist length. Let us as- 
sume that the electric field within the cylinder is equal 
in absolute value to that in the incident wave, i.e. that 
in the absence of the cylinder (the field penetrates the 
conducting medium to skin-depth). In general, this as- 
sumption is wrong and cannot be defended. However, 
in some important cases it does hold true and allows to 
obtain a simple estimate and understand the scaling of 
the problem. Let us also assume that within the beam 
waist the electric field in the incident wave is equal in 
absolute value to that in the center of the beam waist E 
(this assumption cannot lead to additional overvaluation 
by more than e times). The beam power equals 



W = S 1 \P\ = nrj. 1 -.^ 



(103) 



Here Si is the area of the transverse section of the beam 
waist, and |P| is the averaged Poynting vector in the 
beam waist. The power absorbed in the cylinder equals 



2 c 1 2 2-kuj ■ 4a 2 



Here V = Tra 2 li is the volume where absorption takes 
place, 



a one-wave approximation. The idea is as follows. Any 
incident beam can be expanded in cylindrical waves (see 
Eqs. H78J1 and l|79jl ). For each cylindrical wave there is a 
relatively simple exact solution of the problem of diffrac- 
tion on an infinite cylinder. Let us choose a typical cylin- 
drical wave from the expansion. Its field on the axis will 
have a structure similar to that of the incident beam. We 
may then assume that for the incident beam within the 
beam waist length the refracted wave coincides with that 
for the chosen cylindrical wave (the latter should be nor- 
malized in such a way that the magnitudes of the fields 
in the incident beam and in this wave are approximately 
equal). So let us assume that the actual incident, re- 
flected, and refracted fields may be approximated within 
the beam waist by Eqs. (|81|) . l|82|) . and JSSJl, correspond- 
ingly. We further assume that subscript n — 1 (otherwise 
there is no efficient absorption in the longitudinal geome- 
try, at least if Ai <C 1, which is true for cylindrical waves 
that have a field structure similar to that in the waist of 
Gaussian beams), and this subscript will be omitted. Let 
us choose 



M & —> 



(106) 



so that the width of the Gaussian beam waist is approxi- 
mately equal to the width of the main maximum of func- 
tion J (Aip). We are interested in analytical conditions 
of efficient absorption of broad beams in thin cylinders, 
so we make the following assumptions: 



Pi < 1, 



PI 



> 



1 J'(P2) 



P2 J(P2 



pi 



> 



£2 1 J'(P2) 



£l P2 J(P2) 



\P2\ »Pl, 



(107) 

(108) 

(109) 
(110) 



(105) 



is the beam waist length for the Gaussian beam 
(Ref. 14]). Thus, the absorption efficiency equals 
unity to the order of magnitude. One can see that if 
the beam waist width increases twofold, the electric field 
in the cylinder falls fourfold (if beam power does not 
change), but the beam waist length increases fourfold, 
so the absorption efficiency does not change. Alterna- 
tively, if we use a cylinder with the radius that is half 
as large, the optimum conductivity increases fourfold, so 
in this case the absorption efficiency does not change, 
either. However, this estimate calculation is not always 
correct, so let us describe a semiquantitative method of 
calculation of the absorption efficiency. This is actually 



Ai < 1. 

Eq. (|111|) implies 

7 wl. 

From here on we assume that 

Ml = M2 = £i = 1 



(111) 
(112) 
(113) 



(the ambient space has the properties of vacuum, and 
the material of the cylinder is non-magnetic) and omit 
the subscript in S2- Then, taking into account that 



2i 

7TX ' 



(114) 
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(115) 



for x <C 1, we obtain: 



a2 ~ ae 



1 



-2i 



1 



-2 1 



1 



D V Pi 

i 



-is + 



ie(a + i0), 
2 1 



+lj pp(p 2 ) Pl D { - £) = pp(p 2 )piD 



b 2 « a 



1 



J (pz)it P 2 2 {-^-\DPl 



Tl 



+ 



+0 



1 



-2i 



J (P2) Trpl^-^ D \ Pi) Pp(P2)piD 



1 \ i 
a- 



pp(p 2 )piD p 2 2 J(p 2 )piD 



(a + 10). (116) 



Let us estimate power in the incident beam. If complex 
amplitudes of electric and magnetic fields E am and H am 
in the center of the beam waist (p = 0, z = 0) are the 
same as those for the linear combination of cylindrical 
waves of Eq. JHU, then we obtain using Eqs. (f?%|) . fT^jl : 

£„ m «(a+i/3)y{i,-l,0}, 

ff„ w «(a + */3)y{l,» > Q} (117) 
(we take into account that 



lim 



p^O p 



2 ' 
■/i(0) = ~ 



(118) 



use Eqs. (|lll|l . (|112|l . and omit factor exp(my)). The z- 
component of the Poynting vector in the center of the 
incident beam averaged over a period equals 



1 c 

2 4^ R ' e [- EQm X H am\* 



1 C 



2 ^ R&X-EampHamtf -^avatpH amp ) — 

= ~ Re((a + i/3)(a' - i/T)^ • 2 = 9^-, (119) 
2 47T 4 1d7t 

where 9 — aa* + 13(3* + ia*f3 — ia.0*. We may assume 
that the effective area of the transverse section of the 
Gaussian beam waist equals 



/ 



P 

dxdyexp( -) = nrf, 



(120) 



where r± is the radius of the beam waist, then power in 
the incident beam equals 



W = e^r\=e C ^ = 6 C J^. (121) 
lbn 16 16ar 



Let us estimate the energy flow through the lateral sur- 
face of the cylinder. The p-component of the Poynting 
vector at z = 0, p = a averaged over a period equals 



l_c_ 

24^ 



Re[£#>, <p, 0) x H* r *(a, <p, 0)] p . (122) 



From Eqs. (|78|l . (|79(l . and (|83|) we find for complex am- 
plitudes of the fields in the refracted beam: 

H£ r {a,<p,a)=a 2 H +b2H( ) , 
E$ r (a,tp,0)=a 2 E +b 2 E' , (123) 

where 



H = -{J(p 2 ),ip 2 J'(p 2 ),0}, 
a 

E = -< J (p2), — J{p 2 ),—J{p 2 ) 

aye e ea 

R' = - ihP2J'(p2), -lJ( P 2), ^J(P2)) , 

K = ~ W(pa). -itoJ'fa), 0} (124) 



(we omit factor exp(iny?)). Therefore, in view of 

Eq. nm 



E$=E$(a,<p,0) 



%<p dip 

(a + i0)—j^ — - fie f ] J +{-ip 2 J') 



ap 2 Jp\D 



{a + i0) 



-A-iKJ+p 2 J'), 



ap\Jp\D ' 



{a + i0)— j- — —ie—J = 
ap 2 Jp\D ea 

= (a + if3)-^— i ^J, 
ap 2 Jp\D a 

H% = H%(a,<p,0) = 



(a + i0) 



ap 2 Jp\D 



(a + i0) 



ap 2 Jp\D 



(ie ip 2 J' + (—J)) = 
(-l)(J + ep 2 J'), 



H r J z r = H^{a^.,Q) = 

= (a + i0 1 glj, (125) 
ap 2 Jp\U a 
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where J = Ji(p2), J' = J[(j>2)- 

Re[^(a,^0)xJf^*(a,^0)] p = 

= (a + i/3)(a* -0*)- 



'■\P2W\ 2 P\\D\ 2 



xRe i(J* + p* 2 J'*)^J - J(-1)(J* + £*j>2 J'*) 
a a 



a 3 | P2 | 4 P?|£| 2 



7'* 

Im [pi ( — -e*p* 



(9 



a^pW 



1* 

Imlpf (2 + (e+l)p 2 ^) ) = 



J 



20 



(126) 



more accurate estimates. Let us introduce the following 
designations: 



Then 



h(pi) = 

fa = faip-2) 



1 H'{ Pl ) 
Pi H( Pl ) ' 
, 1 J' fa) 

P2 J{P2) ' 



D 



Pi 



pi 



l>2 



_ -Cfi(Pi) - h{P2))(h{Pi) ~ eh(P2)). 
Let us introduce an additional assumption: 

|ln(pi)| »!• 



(130) 



(131) 



(132) 



To estimate the energy flow through the surface of the 
cylinder within the beam waist, we should multiply the 
Poynting vector of Eq. I|121(l by the area of the lateral 
surface of the cylinder: 



2irali 



(127) 



We neglect the energy flow through the sections of the 
cylinder at z — ±y, as in view of Eqs. H105|) . l|106|) 
it is much less than the power in the incident beam 
(Eq. (|121jn . The power absorbed in the cylinder W a 
equals the energy flow through the surface of the cylin- 
der with an opposite sign (as positive p-component of the 
Poynting vector corresponds to outward energy flow): 



o , 1 c 29 T ( 1 e + l J' 

W a = -2-Kah ^r, — r^Im H 

1 2 4tt a 3 p 2 \D\ 2 \pl 2 p 2 J 



Trk 
' 2tt 



s2 C 20 T 

Aa 2 p{\D\ 2 



1 e + l J' 



p\ 2 V2J 
6cr 2 t (I e + l J 1 



a 2 p\\D\ 2 ~ 



P\ 2 P2J 



.(128) 



We used Eqs. {TUB) . fT2li|) . From Eqs. (JUJl.CJl we ob- 
tain an expression for the part of the power of the incident 
beam that is absorbed in the cylinder: 



W a _ 

~W ~ ~ p\\D\ 2 



16 I e + l J 1 
Im ( — + — 

P2 2 P2J 



(129) 



Let us estimate determinant D (Eq. (|97|0 . We use the 
assumptions of Eqs. (|107I113[I . but sometimes we need 



This assumption may seem unrealistic, as it means that 
the radius of the cylinder is several orders of magnitude 
less than the radius of the beam waist. However, this as- 
sumption allows to derive asymptotic formulae that give 
very useful guidelines for more realistic cases. It seems 
worthwhile to analyze the case of extremely thin cylin- 
ders, as they are more difficult to heat. Thus, let us use 
the following expansions for the Hankel functions with 
small arguments: 



7T V 2 X 



7T \2 



(133) 



so 



1 3 ln Pi + ^ 1 l + ip?ln 



Pi mpi 



1 



-^(l+p?lnpi).(134) 
Pi 



Using Eq. (|SU|) . we obtain: 



■^ = --l'-g 



2 2 22 

P2 -, Pi -, P2 - Pi 

a a a 



(135) 



Therefore, 
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D 



2 

pi 



2 

Pi 



v\-v\ 



1 



2 2 

pip 2 



i 



i 



2 
P5 



1 P| 
-9 - 2/3 — fo—z 



\sxp x + f 2 (p2 
1 2 



^ + lnpi+/2(p 2 ))/ 2 (p 2 )^-/i 
Pi / « 



a 2 p 2 



hi'. 



vl 



pi 



1 



a?p\ 



1 



(136) 



where 



where 



/ 3 = lnpi + / 2 (p2)- 



(137) 



We further assume that the absolute value of the term 
with in the last line of Eq. I|13tjfl is greater than or 

of the same order of magnitude as the absolute values 
of the two other terms. As D enters the denominator 
of Eq. (|128|) . we may now neglect the two lesser terms: 
for any domain of parameters providing high absorption 
efficiency that we may find in this way, the assumption 
will ensure high absorption efficiency after we take into 
account the two neglected terms. On the whole, our as- 
sumptions may be loosely defined in the following way: 
pi (which is proportional to the ratio of the radii of the 
cylinder and the beam waist) is "much less than any other 
parameter of the problem" . Thus, we obtain 
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and 
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-2 + h 1 

P2 



A 
2a 2 



1 



k = h - m Pi, h = --r = 
h 



Therefore, 



Pi 



h 1 



2a 2 



-/ 2 1 



2a 2 
/a 



P 2 



2 

P2 



1 



1 

~ 2a 2 _ 
2a 2 

2 

2a 2 



(140) 



(141) 



On the other hand, 
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Thus, from Eqs. (|139() . (|140|l . (|140|> . we obtain the asymp- 
totic formula for absorption efficiency: 
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(pf + 2lP 



In pi 



(143) 



We see that, all other parameters being equal, the ab- 
sorption efficiency depends very weakly on p±, which is 
proportional to the ratio of the radii of the cylinder and 
the beam waist. 



V. EFFICIENT HEATING IN THE 
TRANSVERSE GEOMETRY 

Let us now consider the transverse geometry, where 
a cylindrical electromagnetic wave converges on a thin 
conducting cylinder. The general solution of Section ITTT1 
may be used in this case with (3 = 0, n = 0, 7 = 0, but 
this relatively simple case is worth direct consideration, 
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as it is easier to verify, and it is more natural to choose 
the incident TM wave using the Hankel function of the 
second kind rather than the Bessel function, as the inci- 
dent wave propagates in the radial direction (towards the 
axis) rather than along the axis. Converging TE waves 
are not considered here, as they do not provide efficient 
heating in the relevant parameter domain. 



IT 



wpinc 



{o,^ 2) '( P ),o}, 

= {0,0,i4 2) (p)}. (144) 



(2) 

Here H n is the Hankel function of the second kind. We 
assume here that p,\ = p, 2 = 1, E\ = 1, e 2 = £,fco = 1- 
Again, the actual electric and magnetic fields are equal 
to the real parts of the products of the above expressions 
by the time-dependent factor exp(— iuit). A diverging 
reflected cylindrical wave propagates from the cylinder 
outwards: 

H^ 1 =a 1 {0,z^ 1) '(p),0}, 

oi{0,0,^ 1) (p)}. (145) 
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Inside the cylinder, the refracted wave has the following 
form: 

n r/r = a 2 J (Vip), 
H* r = a 2 {o,iViJo(Vtp),o}, 
E^ T = a 2 {0,0, Jo(v^p)} • 



(146) 



The tangential components of the electric and magnetic 
fields must be continuous at the surface of the cylinder, 
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iH^ (a) + iaxH^p (a) = ia 2y /eJ (y/ea), 

H { 2) (a) + aiH^ia) = a 2 J Q {^a). (147) 
Hence, if we denote \fea = p, then 
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Using the formula for the Wronskian of Hankel functions, 
we obtain: 
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The Poynting vector averaged over a period can be ex- 
pressed in terms of complex amplitudes of electrical and 
magnetic fields as |j^Re[-E x H*] (the asterisk denotes 
complex conjugation). The heating efficiency rj (the share 
of the power in the incident wave that is absorbed in the 
cylinder) equals the ratio of the radial components of en- 
ergy flows in the refracted wave (Eq. I|146jl ) and in the 
incident wave (Eq. (|144fl ) at the surface of the cylinder. 
These components do not depend on the azimuthal angle 
if, so 

Re[E r f r x H r f r *] p _ Re{-E r J r H r J r *) 
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(151) 



We are interested in the parameter domain where the 
conductivity of the cylinder is high and its radius is much 
less than the wavelength. So let us assume now that 
e — e' + ie" , where s' and e" are real, e" 3> |e'|, and 



a < 1, \p\ < 1, I lna| > 1. 



(152) 



Then Up) » 1, j' Q (p) « -§, H^(a) « iflna, 
• ff o 1) ( a ) ~ ^1 so E< 4- (SB) becomes: 
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On the other hand, 
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Thus, 
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At 
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(150) 



If 
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(157) 
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where ~ means " is of the same order of magnitude as" , 
then 

(158) 



\na\ £ " a2 l lna l 
1 2 



2 | lna| 



Unal 



It is evident that the heating efficiency decreases ex- 
tremely slowly when the radius of the cylinder decreases 
by many orders of magnitude, if in the same time the 
conductivity of the cylinder increases to satisfy the opti- 
mality condition (Eq. I|157j l). In this case the efficiency 
typically equals 20 ~ 40%. Moreover, the condition of 
high efficiency is not resonant, as the width of the maxi- 
mum corresponds to variation of conductivity by an order 
of magnitude (at the level of over 50% of the maximum). 

It should also be noted that Eq. (|157Jl is compatible 
with the approximations (Eq. (|152H ). 



VI. CONCLUSIONS 



greater than those of the cylinder. The relevant condi- 
tions are not too rigid. 

This is the first part of the work. It contains the fol- 
lowing results: 

Exact solutions of the free Maxwell equations are de- 
scribed that are approximated by Gaussian beams with 
great accuracy when the beam waist radius is much 
greater than the wavelength, and their asymptotic prop- 
erties are proved. It is also proven that these solutions 
are normalizable. These results may be useful for many 
applications. 

The qualitative, semiquantitative, and quantitative ap- 
proaches are developed for calculation of efficiency of 
heating of cylindrical targets with broad electromagnetic 
beams. The quantitative approach includes rigorous so- 
lution of the Maxwell equations. 

Asymptotic formulae for heating efficiency in the lon- 
gitudinal and transverse geometry are derived. 



A rigorous analysis of a diffraction problem demon- 
strates that it is possible to achieve efficient heating of 
cylindrical targets by electromagnetic beams with trans- 
verse dimensions that are several orders of magnitude 
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